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ABSTRACT 

We present an analysis of deep XMM-Newton and Chandra observations of the z = 
1.05 galaxy cluster XLSSJ022403.9-041328 (hereafter XLSSC 029), detected in the 
XMM-Newton large scale structure survey. Density and temperature profiles of the 
X-ray emitting gas were used to perform a hydrostatic mass analysis of the system. 
This allowed us to measure the total mass and gas fraction in the cluster and define 
overdensity radii -R500 and i?25oo- The global properties of XLSSC 029 were measured 
within these radii and compared with those of the local population. The gas mass 
fraction was found to be consistent with local clusters. The mean metal abundance 
was O.lSto'isZQ, with the cluster core regions excluded, consistent with the predicted 
and observed evolution. The properties of XLSSC 029 were then used to investigate 
the position of the cluster on the M — kT, Yx — M, and Lx — M scaling relations. In 
all cases the observed properties of XLSSC 029 agreed well with the simple self-similar 
evolution of the scaling relations. This is the first test of the evolution of these relations 
at z > 1 and supports the use of the scaling relations in cosmological studies with 
distant galaxy clusters. 

Key words: cosmology: observations ~ galaxies: clusters: general - galaxies: high- 
redshift - intergalactic medium - X-rays: galaxies 



1 INTRODUCTION 

The mass function of galaxy clusters is predicted by theo- 
retical models of the formation of structure from the den- 
sity fluctuations in the early universe. The shape and evo- 
lution of the predicted mass function are strongly depen- 
dent on the details of the input models, and so compari- 
son of predicted and observed mass functions at difl'erent 
redshifts can place tight constraints on the values of inter- 
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esting cosmological parameters (e.g. [ Reiprich fc Bohringeij 
l2002l : ivikhlinin et al.|[20o3 : lHenrvll2004l '). However, in order 
to perform these tests, a key observational challenge must be 
met; galaxy cluster masses must be determined from their 
observable properties. This is an area of active research, and 
many approaches have been investigated using different ob- 
servables, including the properties of the hot ionised intra- 
cluster m edium (IC M) as measur ed from its X-ray emission 
ISarazinl 198S:,Rosati et aLll2002l ') or the Sunyaev-Zel'dovich 
effect I Birkinshawl 19991 ): the distribution of properties of 
member galaxies such as richness and velocity dispersions 
and the strong and weak gravitational lensing effects of 
cluster gravitational potential s on background galaxies (e.g. 
ISmith et al1l200ll : lDahlell2006l) . Numerical simulations of in- 
dividual galaxy clusters and large cosmological volumes are 
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also an invaluable tool for studying clust ers and testing these 
different observat ional techniques (e.g. lEvrard et al.l ll99"^ : 
iNagai et al.ll200'iD . 

X-ray observations of clusters are useful because the 
ICM is extremely X-ray luminous, allowing the detection 
and study of galaxy clusters out to 2;>1. Under the as- 
sumption that the ICM is in hydrostatic equilibrium with 
the cluster's gravitational potential (a reasona ble assump- 
tion for clusters that are not actively merging; IPoole et al.l 
I2OO6I ) then radial profiles of the gas density and tempera- 
ture measured from X-ray data can be used to determine 
the total cluster mass within some radius. The radius used 
is typically chosen to enclose an overdensity A with respect 
to the critical density at the cluster's redshift, with A = 200 
approximating to the virial radius, and A = 500 the maxi- 
mum radius detectable in typical observations. The hydro- 
static mass estimation described above demands high qual- 
ity X-ray data with which to measure the required profiles. 
Such data are becoming more commonplace in the Chan- 
dra and XMM-Newton era, but are still far from the norm. 
In general, and particularly for distant clusters (which pro- 
vide the most information for measuring cosmological pa- 
rameters). X-ray data permit the measurement of simple 
global properties such as a single gas temperature and lu- 
minosity. Such properties can still be useful for estimating 
cluster masses, as power law scalin g relations exist between 
globa l properties and cluster m ass jKaisedllQsd : Markcvitch 
1 19981 : iFinoguenov et ahl I2OOII ) . These scaling relations are 
predicted by simple self-similar models of clusters. While 
the observed relations can differ from the self-similar predic- 
tions, the slope, normalisation and evolution of the relations 
can be measured using high quality data, and then applied 
to give mass estimates for poorer quality data. These tech- 
niques have allowed studies using the temperature and lu- 
minosity functions of galaxy clusters as cosmological probes, 
with scaling relations used to convert between the observ- 
ables and masses for comparison with theoretical mass func- 
tions. 

The scaling relations are, however, an imperfect tool 
for converting observables to masses. Radiative cooling and 
cluster mergers cause significant scatter in the relations be- 
tween X-ray luminosity (I/x), temperature (kT) and mass, 
limiting the accuracy of derived masses. Much of the scat- 
ter derives from the dense core regions of clusters where 
radiative cooling can be extremely efficient, and the effects 
of mergers are most pronounced. Excluding t hese regions 
can significantly reduce the resulting sc atter (|Markevitchl 
19981: lO'Hara et al.ll2006l : iMaughanI |2007| ) . Recent work by 



Kravtsov et al.l 1 20061 ) has shown that the parameter Yx , the 



product of gas temperature and gas mass has an extremely 
low scatter with total mass, providing the potential for reli- 
able cluster mass estimates from simple observables. 

The mass-observable scaling relations are reasonably 
well-calibra ted for relaxed clu s ters in the local un iverse 
{z < 0.15) l|Arnaud et al.ll2005l : Ivikhlinin et al.ll2006l ), but 
at higher redshifts, the relations are less well measured. This 
is due to the long observations required to obtain sufficiently 
deep X-ray data to allow hydrostatic mass estimates against 
which to plot the simple observables. Recent work, based on 
small samples of clusters with the best available data has 
enabled the measurement of th e M — kT relation at z < 0.7 
IIKotov fc Vikhlininll2005l ,[ 20061 ). The relations were found to 



have the same slope as their local counterparts, with the evo- 
lution of the normalisation well described by the s elf-similar 
model. Furthermore, using some of the same data, iMaughanI 
(2OQ7I ) found that the Yx — M relation at 2; ~ 0.6 is consis- 
tent with self-similar evolution of the local relation. Beyond 
z ~ 0.7, X-ray hydrostatic mass estimates have been made 
for only one cluster CLJ1226. 9-1-3332, although that mass 
estimate may be bias ed by merging activity in the cluster 
('Ma ughan et al.ll2007l ). The evolution of the mass-observable 
relations to such high redshifts have thus not been well stud- 
ied, although more a nd more clusters are being detected at 
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these dista nces (e.g. iRosati e t ai?'200- 
[Stanford et al. 200^; iBremer ct al. 20oi; 
providing a powerful resource for cosmological studies, given 
well-measured scaling relations. In this paper we make an 
important step forward in this endeavor by measuring the 
first X-ray hydrostatic mass for a cluster at z > 1 and inves- 
tigating its position on the various mass-observable scaling 
relations. 

Galaxy cluster XLSSJ022403. 9-041328 (hereafter 
XLSSC 029) was detected as an extended X-ray source in 
the XM M-Newton Large Scale Structu re (XM M-LSS) sur- 
vcy ( Andreon et al. I I2005I : n ierre et al.i ,2006: Pacaud et al 



2007jV The XMM-LSS survey is described in iPierre et al" 



(2004 I) , and som e results from the survey ar e presented in 
FValtcha nw et al.l (|2004); IWiUis et al.l ^M^; (Pierre et all 
(200fj); IPacaud et al.l (|2007l ). NIR imaging confirmed an 
overdensity of faint galaxies coincident with the X-ray 
source XLSSC 029, and spectroscopic follow-up confirmed 
the redshift of the cluster to be z = 1.05. The cluster 
has since been the target of deep Chandra ACIS-S and 
XMM-Newton follow-up observations, which are the subject 
of this paper. In the following sections we describe the 
reduction of the X-ray data, present the results of our mass 
analysis, and investigate the location of the cluster on the 
different mass-observable relations. A ACDM cosmology of 
Ho = 70 km s"^ Mpc"^ ) = 1, and Qm = 0.3 {Ha = 0.7) 
is adopted throughout and all errors are quoted at the 
68% level. In this cosmology, at z = 1.05, 1 arcsecond 
corresponds to 8.1 kpc. 



2 X-RAY DATA REDUCTION 

The reduction an d analysis methods use d were the same as 
those presented in iMaughan et~al] (|2007 l) for both the Chan- 
dra and XMM-Newton data, but the most important details 
are repeated here. The standard data reduction procedures 
were followed, using the most recent calibration products 
available as of March 2007. Lightcurves were produced for 
the observations, and were cleaned to remove periods of high 
background. The Chandra observation of XLSSC 029 was 
taken in six consecutive exposures with cleaned exposure 
times of between 10 ks and 30 ks. The total Chandra expo- 
sure time was 127 ks, yielding a total of ~ 1300 net source 
counts in the 0.3 — 5 keV band. The XMM-Newton observa- 
tions had a cleaned exposure time of 70 ks (pn) and 88 ks 
(MOS) giving a total of ~ 3400 net source counts in the 
0.3 — 5 keV band. The observations used are summarised in 
Table [U 



© 0000 RAS, MNRAS 000, 000-000 



XMM-Newton and Chandra observations of XLSSJ022403. 9-041328 3 



Date 


Detector 


Obs ID 


Total Exposure (ks) 


Good Time (ks) 


iUUo-Ui-Ui 


XMM-Newton pn 


OZiU4yOiOi 


86 


70 


2005-01-01 


XMM-Newton MOSl 


0210490101 


106 


87 


2005-01-01 


XMM-Newton MOS2 


0210490101 


106 


88 


2005-09-13 


Chandra ACIS-S 


6390 


12 


10 


2005-10-12 


Chandra ACIS-S 


6394 


18 


18 


2005-10-12 


Chandra ACIS-S 


7182 


23 


23 


2005-10-14 


Chandra ACIS-S 


7184 


23 


23 


2005-10-15 


Chandra ACIS-S 


7183 


20 


20 


2005-11-21 


Chandra ACIS-S 


7185 


33 


33 



Table 1. Summary of the XMM-Newton and Chandra observations of XLSSC 029. 



2.1 Background preparation 

For the analysis of both Chandra and XMM-Newton data, 
background estimates were derived from blank-sky datasets, 
using the same detector regions as the source emission being 
considered. Important differences exist between the blank- 
sky data and the observations of a given source, due pri- 
marily to differences in the levels of the particle-induced 
background components (which can vary significantly with 
time and dominate at energies > 2 keV) and differences in 
the soft Galactic foreground emission (which can vary signif- 
icantly with position on the sky and dominates at energies 
< 2 keV) . These differences must be corrected for when us- 
ing the blank-sky data. 

In the case of XLSSC 029, due to the target's high red- 
shift, the source observations contained large detector re- 
gions that were free from source emission. Spectra were ex- 
tracted from these "local" regions and compared with the 
blank-sky spectra from the same detector regions. 

For the XMM-N ewton data, the blank-sky files of 
ICarter fc ReadI l|2007l ) were used. Periods of high particle 
background were removed by filtering lightcurves produced 
in the 10 - 12 keV (MOS) and 12 - 14 keV (pn) energy 
bands with the same count rate limits that were applied 
to the source data when they were cleaned. The blank- 
sky spectra were then normalised to match the local back- 
ground count rates in these high energy bands and com- 
pared with the local background spectra. Initially, a very 
poor agreement was found, and additional flare filtering at 
lower energies (2 — 5 keV, using iterative 3cr clipping of the 
lightcurve) was required for the blank-sky files before there 
was a good agreement between the local and blank-sky back- 
ground spectra. This is illustrated for the pn camera in 
Fig. [T] and very similar results were found for both MOS 
cameras. After this additional filtering, the ratios of the en- 
tire field local to blank-sky count rates in the high energy 
bands were 1.01, 1.10 and 1.14 for pn, MOSl and MOS2 re- 
spectively, indicating that the particle background levels are 
similar in the target and blank-sky fields. After scaling by 
these factors, the XMM-Newton local and blank-sky spectra 
agreed well except below ~ 2 keV where there was a signif- 
icant decrement in all local background spectra compared 
to the blank-sky spectra. This is due to significantly below 
average soft Galactic foreground emission in the direction of 
XLSSC 029 and is corrected for in subsequent spectral and 
image analysis. 

For the Chandra observations, the standard blank-sky 



background files were usecQ. The blank-sky spectra were 
then normalised to the local background spectra in the 
9.5 — 12 keV band to account for variations in the parti- 
cle background and compared. A good agreement was found 
in all of the Chandra observations, with no additional filter- 
ing required. Consistent with the XMM-Newton background 
spectra, there was a significant decrement at < 2 keV in all 
local background spectra due to the low soft Galactic fore- 
ground emission in the direction of XLSSC 029. 



3 IMAGE ANALYSIS 

Images were produced for each dataset in the 0.3 — 2.0 keV 
(observed frame) energy band, and the Chandra and XMM- 
Newton data were co-added to make combined images for 
each satellite. These images were then vignetting corrected 
(using an exposure m ap generated for 1.5 keV photons) and 
adaptively smoothed l|Ebeling et al.l[200 6^). Contours of the 
smoothed Chandra data are plotted in Fig [21 overlaid on 
an image produced from a 2400 s F0RS2 VLT I band ex- 
posure taken in 2005 August (proposal ID 075.A-0175, PI 
Andreon) with 0.7 arcsecond seeing. The contour levels were 
defined so that the emission bounded by adjacent contours 
was detected at a significance of 3a above that boun ded 
by the surrounding contour (see iMaughan et al. I l2007l . for 
details). The X-ray morphology is slightly elliptical, but ap- 
pears reasonably relaxed, with no indication of substructure 
in the images. The X-ray morphology in the XMM-Newton 
image was consistent, with stronger contamination from 
point sources due to the larger point spread function (PSF). 
The centroid of the X-ray emission (determined from the 
Chandra data alone) is located at a[2000.0] = 02''24™4.0° 
(5[2000.0] = -04°13'28.9". 

Wavelet de composition was used to detect point sources 
in the images (|Vikhlinin et al.l 1 19981 ) , and detected point 
sources were excluded from all analysis of the extended emis- 
sion. There was no detection of point-like emission from the 
central galaxy. The X-ray point sources are apparent in Fig. 
[2] which also shows the position of a faint source detected 
by XMM-Newton but not by Chandra. Finally, a low fre- 
quency radio source was detected at 325 MHz (J0224. 0-0413 
at a[2000.0] = 02''24'"3.82" (5 [2000.0] = -04°13'3 8.8"), close 
to the core of XLSSC 029 bv lCohen etlo] l|2003h . No coun- 



|http : / /cxc . harvard . edu/contrlb/maxlm/aclsbg/ | 
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Figure 1. Comparison of the local (black points) and blank-sky (grey points) background spectra measured with the XMM-Newton pn 
camera. The two plots show the effect of filtering the blank-sky data for background flares in the high-energy band alone (left plot) and 
of additional filtering in the 2 — 5 keV band (right plot). 
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Figure 2. Contours of the adaptively smoothed image produced from the Cha ndra data ar e overlaid on a VLT I-band image of XLSSC 
029. The dotted circle shows the position of the 325 MHz radio source detected byl Cohen et al. (2003) with no apparent X-ray counterpart. 
The cross near the bottom right indicates the location of a faint X-ray point source detected in the XMM-Newton data that was not 
detected by Chandra. 



terpart for this source was detected in either the Chandra 
or XMM-Newton images. 

Using the Chandra image alone, due to the smaller 
PSF, the morphology of the X-ray emission was analysed 
quantitatively. A two dimensional /3-model was fit to the 
X-ray image (including instrumental effects) and the el- 



lipticity of the best-fitting model was e — 0.28 ± 0.05. 
The centroid shift (the standard deviation of the separa- 
tion between X-ray peak and centroid) was measured to be 
(w) — (0.012 ± 0.001)-R5oo {R500 is the radius correspond- 
ing to A = 500 as computed in i|4.ip . These centroid shift 
and ellipticity values are close to the median values found 
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Figure 3. XMM-Newton surface brightness profiles of XLSSC 
029 and tlie blank-sky background. The profiles were produced in 
the 0.3 — 2 keV band and the background profile was normalised 
to match the source profile in regions free of source emission. 
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Figure 4. The emission measure profile of XLSSC 029 measured 
from the Chandra and XMM-Newton data are plotted along with 
the best fitting model. The lines show the same 3D model pro- 
jected along the line of sight, and (in the case of XMM-Newton) 
convolved with the telescope PSF. 



for a sample of 115 gala xy clusters of a range of redshifts 
observed with Chandra (jMaueha 3 12OO7I). If only the more 
distant clusters {z > 0.5) from the iMaughanI (|2007l ) sam- 
ple are considered, then the centroid shift of XLSSC 029 is 
lower than 21/32 of those clusters. XLSSC 029 can thus be 
considered a relatively relaxed example of a distant cluster. 

In order to determine the mass of gas in the cluster, 
a surface brightness profile was extracted from each image 
and converted into an emission measure profile (as described 
in [Maughan et al. 2007.) . This conversion assumed a metal 
abundance of O.SZ©, and a temperature that varied with 
projected radius accordin g to th e mean temperature profile 
found by IVikhlinin et~al] l|2006l ). normahsed to the global 
temperature and i?5oo measured for XLSSC 029 ( ti4.2|l . The 
conversion is only weakly dependent on temperature so the 
choice of temperature profile does not significantly afi'ect our 
results. 

The background level in each bin was estimated from 
images produced from blank-sky background files that were 
normalised to match the count rates in the imaging energy 
band in regions free of source emission. This was an iter- 
ative process, with the extent of the cluster emission be- 
ing determined from the profiles and then the background 
images renormalised excluding this emission, with the pro- 
cess repeated until convergence. The background normalisa- 
tion is illustrated in Fig. [31 which shows the radial surface 
brightness profiles of the source and background (after vi- 
gnetting correction) for the combined XMM-Newton images. 
Vignetting correction of the profiles was performed by divid- 
ing the counts in each bin by the sum of the exposure map 
counts in that bin. Note the upward curve of the profiles at 
large radii is due to the boosting of non-vignetted particle 
events by the vignetting correction. This background com- 
ponent is correctly subtracted, however, because the particle 
events are boosted in the same way for the source and back- 
ground images, and level of particle events is similar in the 
source and background datasets (^STTJ. The good agreement 
of the source and background profiles outside of the cluster 
emission indicate that the blank-sky background files pro- 
vide an accurate background level for the emissivity profiles. 

A three-dimensional (3D) model emission measure pro- 



file was then projected along the line of sight and fit si- 
multaneously to the observed Chandra and XMM-Newton 
profiles. For the XMM-Newton data, the projected profile 
was convolved with a model of the XMM-Newton PSF be- 
fore the fit statistic was computed. The PSF model used was 
a King functiorQ with a core radius of 5.4" and a slope of 
1.5. These parameters are appropriate for 1.5 keV on-axis 
photons and are the mean of the MOS and pn values. 

The emission measure mod el used was the modified l3- 
profile of lVikhUnin et all (|2006l ): 

-^"^ - -0(^^^^J;|5r^(l + rVrJ)-^/-. (1) 

The best-fitting model is plotted in Fig.|4]along with the ob- 
served profiles. The model was found to be a good fit to the 
combined data, with x^/^ = 45.8/43 demonstrating good 
agreement between the Chandra and XMM-Newton data. 
The best model parameters were no = 1.31 x 10~^ cm~^ 
rc = 55.7 kpc, a = 0.0, /? = 0.46, = 872 kpc, and 
e — 5.0 with 7 fixed at 3. Errors were not computed 
on individual parameters as they are highly correlated, 
and model fits to Monte Carlo randomisations of the data 
were used to determine the uncertainties on the gas mass 
and other derived properties. The gas density is related 
to the emission measure by pg — 1.252mp{npney^^ , as- 
s umin g ja, cosmic plasma with helium abundances given by 
Aiidersji Grevesso (1989). 

While the morphology of XLSSC 029 is elliptical, our 
analysis assumes spherical symmetry. However, this assump- 
tion leads to negligible errors ( < 3%) in the total and 
gas mass compared to that derived using a triaxial model 
(Piffarctti ct al. 2003). To allow comparisons with other 
clusters where such detailed modeling of the gas distribution 
is not possible, the Chandra s urface brightness profile was 
fit w ith a standard /3-model (jCavaliere fc Fusco-Femiand 
Il976h . The best-fitting model had a core radius of rc = 
84 ± 5 kpc and a slope of /3 = 0.60 ± 0.02. 

^ |http : //xmm. vllspa. esa. es/docs/documents/CAL-TN-QOlS .pdf] 
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4 SPECTRAL ANALYSIS 

The general procedure followed for the analysis of the X- 
ray spectra of XLSSC 029 was to extract spectra in a se- 
ries of annular bins, and fit these to create a temperature 
profile. The temperature profile was then used to deter- 
mine the mass profile of the cluster and define R500. Spec- 
tra were then extracted from within this radius and fit to 
give the global spectral properties of the syste m. All spec- 
tra w ere modeled by an absorbed APEC model (|Smith et al.l 
I2OO j ). with the absor bing column fixed at the Galactic value 
(2.58 X 10^° cm"2; iDickev fc Lockmaniil99Ql . The model 
temperature, metal abundance and normalisation were free 
parameters in the fits, and the source spectra were grouped 
to contain > 30 counts per bin and fit using minimisa- 
tion. For XMM-Newton the pn and MOS spectra were fit 
simultaneously with the same model, and similarly the data 
from the six Chandra observations were fit simultaneously 
for each region considered. 

For both Chandra and XMM-Newton, background spec- 
tra were extracted from the blank-sky background files us- 
ing same detector regions as the source spectra. The back- 
ground spectra were normalised to match the high energy 
count rate in the source datasets as discussed in ^2.11 In or- 
der to correct for the decrement in soft Galactic foreground 
emission towards XLSSC 029 compared to the blank-sky 
data, soft residual spectra were created by subtracting the 
blank-sky spectrum from the local background spectrum in 
the same, source-free detector region. In the case of Chan- 
dra, the residual spectra were fit with a thermal model with 
negative normalisation. This was then included as a fixed 
component (scaled by extraction area) in subsequent fits to 
the source spectra (th is method is described in more detail in 
IVikhhnin et all 120051 ). For XMM-Newton, the soft residuals 
were scaled by extraction area and added onto the back- 
gro und spectra (a, proce ss known as double subtraction; see 
e.g. lArnaud et al.ll200t ). The XMM-Newton Science Analy- 
sis System task evigweight was applied to all of the XMM- 
Newton source and background events lists in order to ac- 
count for differences in the telescope efi'ective area between 
the source region and the region from which the soft resid- 
uals were derived. The entire analysis was also performed 
using a simple local background spectrum extracted from 
the same dataset as the source spectra, and no statistically 
significant difi'erences were found in any of the derived clus- 
ter properties. 

4.1 Temperature and Mass Profiles 

Spectra were extracted from annular regions centred on the 
cluster centroid defined so that the combined background- 
subtracted counts from the XMM-Newton pn and MOS de- 
tectors were > 500 in each region. The same criterion was 
used to define regions for the Chandra data. The spectra 
were fit and the best-fitting temperatures are plotted in Fig. 
O The temperature profile was then fit with the following 
model; 

This is the same model used bv lVikhlinin et al.l ()2006l ) 
to fit high-quality temperature profiles of local relaxed clus- 



ters but the cool core component has been removed from 
the model as the data for XLSSC 029 do not require it. The 
model is thus a broken power law with a transition region. 
The model was further simplified by setting b = c/0.45, 
thereby fixing the width of the tra nsition region to match 
that found in the average fit to the IVikhlinin et al.l (|2006l ) 
clusters. This 3D model was projected along the line of sight 
to predict the temperature in each bin of the observed pro- 
file. The projection weighted the different temperature com- 
ponents using the measured gas density profil e according to 
the algorithm presented by IVikhlininI (|2006l ) . Due to the 
large radial bins used for the Chandra data, and to simplify 
the projection calculation, only the XMM-Newton data were 
used for fitting the temperature profile. 

The XMM-Newton PSF is not negligible compared to 
the size of the radial bins used for the temperature pro- 
file. The effect of the PSF is to redistribute some photons 
that were emitted in one projected annulus into a different 
projected annulus. The magnitude of this effect was com- 
puted by using a background subtracted Chandra image of 
XLSSC 029 to provide the projected photon distribution 
with no PSF redistribution. For each annular region in the 
temperature profile, a sub-image was extracted, containing 
just the photons in that region. This sub-image was then 
convolved with an image of the XMM-Newton PSF (a PSF 
suitable for 1.5 keV photons detected on-axis was used), and 
the number of photons redistributed from that annulus into 
each other annulus was measured. This process gave a re- 
distribution matrix describing the relative contribution of 
each projected annular region to each bin in the final, pro- 
jected, PSF-convolved profile. In the fitting procedure, the 
temperature in each bin in the final m odel profile was cal- 
culated by a second application of the IVikhlininI (|2006l ) al- 
gorithm to combine the contributions from every projected 
bin, weighted by the PSF redistribution factors. 

The best-fitting 3D model, its projection, and the fi- 
nal model including PSF convolution are plotted in Fig. 
[5] The best-fitting model had parameters fcTo = 3.5 keV, 
rt — 586 kpc, a = 0.20, and c = 1.55. The probability distri- 
butions of the parameters were determined by fits to Monte- 
Carlo randomisations of the data, and used to determine 
uncertainties on derived cluster properties. The best-fitting 
3D model was also compared with the Chandra data plot- 
ted in Fig. [5] The projected temperature predicted by the 
3D model for the inner and outer Chandra temperature bins 
are 4.6 keV and 2.8 keV. The prediction for inner bin agrees 
well with the measured value, while that for the outer bin 
is slightly (but not significantly) higher than the observed 
value {2.2tai keV). 

The best-fitting 3D models for the gas density and tem- 
perature profiles were then used to compute the total mass 
profile of XLSSC 029 under the assumption of hydrostatic 
equilibrium, propagating the Monte Carlo uncertainties to 
all derived properties. This profile was then used to com- 
pute the overdensity profile of the system and determine 
R500 = 0.52t^:J° Mpc and i?2500 = 0.23^0.02 Mpc. The de- 
rived mass profile of XLSSC 029 is plotted in Fig. [U] along 
with the gas mass profile. The total mass of the system 
within i?5oo was found to be 1.3lo;3 ^ lO^^M©. The gas 
mass fraction (/gas) profile of XLSSC 029 is plotted in Fig. 
[71 At the gas mass fraction was 0.14^Q gg, and within 
the smaller radius of -R2500, the gas fraction was 0.08to'.o2- 
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Figure 5. Projected temperature profiles measured for XLSSC 
029 using XMM-Newton and Chandra. The curves sliow the best- 
fitting 3D model, its projection, and the final model after PSF 
convolution, with the la errors on the projected model shown by 
the shaded region. The model was fit to the XMM-Newton data 
alone, and the vertical line indicates -R500 ■ 




R (Mpc) 

Figure 6. Profiles of the derived total and gas mass in XLSSC 
029. The Icr errors on the models are shown by the shaded re- 
gion, the data points mark the midpoints of the XMM-Newton 
temperature profile bins, and the vertical line indicates -Rsoo- 



R (Mpc) 

Figure 7. Profiles of the gas mass fraction in XLSSC 029. The 
1(7 errors are shown by the shaded region and the vertical line 
indicates -Rsoo- 



dra data were not able to provide useful constraints on the 
metallicity. 

The bolometric X-ray luminosity of the cluster within 
J?50o was measured using the data from each satellite. The 
XMM-Newton data gave Lx = (3.4 ± 0.4) x 10'*'' erg s"^ 
(Lx = (4.6 ± 0.4) X lO*" erg s"^ with the central 0.15i?5OO 



'■'^ erg s ^ with the central 0.15-R; 
included), while the Chandra data gave Lx = S.O^ne x 

-HO. 9 in44 



10** erg s-^ 



0.6 

(Lx = 4.5io.4 X 10** erg s"' including the 
core). Uncertainties on the luminosities include those on 
the parameters of the best-fit spectral models. The value 
of Lx given above for the XMM-Newton data with the cen- 
tral O.lS-Rsoo excluded includes a correction for the redis- 
tribution of photons into and out of the excluded region by 
the XMM-Newton PSF. The same method used to calcu- 
lated the redistribution factors for the temperature profile 
was used to calculate that the PSF redistribution leads to a 
net increase of 9% in the number of photons detected in the 
(0.15 — l)-R5oo aperture. The measured luminosity was thus 
scaled down by this factor. 

The global propert ies of XLSSC 029 ar e all consistent 
with those measured by|PacaudelaL| l|200it ) using the orig- 
inal 20 ks of survey XMM-Newton data. 



4.2 Global Spectral Properties 

With -R500 determined from the mass analysis, the inte- 
grated spectral properties of XLSSC 029 were then mea- 
sured within that aperture. The central O.lS-Rsoo was ex- 
cluded, as the emission from this region can be affected by 
strong radiative cooling and merger events, leading to en- 
hanc ed scatter from the self-similar cluster scaling relations 
fe.q. lMaughanll2007l '). The XMM-Newton and Chandra data 
were fit separately. The emission-weighted temperature of 
the gas within -R500 was found to be 3.41q'2 keV based on 
the combined XMM-Newton data, and 4.3lg ^ keV from the 
combined Chandra observations. The m etallicity (iron abun- 
danc es relative to the solar values of lAnders fc Grevessd 
1 19891 ) of the gas was measured to be O.lSlgjgZo from the 
XMM-Newton data. While the data did not permit a metal- 
licity profile of the cluster, the value measured with the core 
regions included increased slightly to OAOt.o'lg'ZiQ, perhaps 
indicative of an abundance peak in the core. The Chan- 



5 SCALING RELATION EVOLUTION 

The hydrostatic mass estimate obtained for XLSSC 029 
presents a unique opportunity to investigate the mass- 
observable scaling relations at z; > 1. Recall that the gas 
mass and total mass were measured using a combination 
of the XMM-Newton and Chandra data, although the to- 
tal mass estimate was strongly dependent on the XMM- 
Newton temperature profile. The global spectral proper- 
ties were measured separately for XMM-Newton and Chan- 
dra, so datapoints for each satellite are included in the 
scaling relations that follow. In all cases, the data were 
scaled by an appropriate power of E{z) (where E^{z) = 
Qmil + zf -I- (1 - f^M - ^^a)(1 + zf -I- Ha, describing the 
redshift evolution of the Hubble parameter) to remove the 
predicted self-similar evolution in each relation. The expo- 
nent used for each relation is indicated in the figures below. 
The mass and temperature of XLSSC 029 were com- 
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Vikhlinin et al (2006) data > 

Vikhlinin et al (2006) fi 

Arnaud et al (2005) data ^ 

Arnaud et al (2005) lit ■ 

XLSSC 029 XMM > 

XLSSC 029 Chandra > 



kT (keV) 

Figure 8. XLSSC 02 9 is plotted on the low -reds hift mass tem- 
peratu re relations of IVikhlinin et alJ l|2006l ) and lArnaud et alJ 
i|2005l '). The two XLSSC 029 points reflect the slightly diflFerent 
global temperatures measured with the XMM-Newton and Chan- 
dra data. The central O.lS-Rsoo was excluded from the tempera- 
ture measurements. 



-4--- 



<4 



Vikhlinin et al (2006) data i — • — i 

Kravtsov (2006) relation 

Arnaud et al (2007) data i — ■ — i 

Arnaud et al (2007) fit 

XLSSC 029 XMM i ' 
XLSSC 029 Chandra ^ 



1e+14 

Yx (Msoi keV) 



1e+15 



Figure 9 . XLSSC 029 is plotted on t he low-redshift Yx - M re- 
lations of IVikhlinin et al.l ( 2006) and lArnaud et al.l ||2TOi). The 
two XLSSC 029 points reflect the slightly different global tem- 
peratures measured with the XMM-Newton and Chandra data. 
The central O.lS-Rsoo was excluded from the temperature mea- 
surements. 



pared with the data of IVikhhnin et all l|2006t ). derived from 
high quahty C handra observa t ions o f low-redshift, relaxed 
clusters and of lArnaud et al] (|2005h . derived from a sim- 
ilar XMM-Newton sample. The properties of XLSSC 029 
were measured in a consistent manner to those of these lo- 
cal samples. The two datapoints for XLSSC 029 (for the 
XMM-Newton and Chandra temperatures) are plotted along 
with the local data in Fig. HJ Both XLSSC 029 datapoints 
are consistent with the local relations. Our data are also 
consistent with the mass- temperature relation measured by 
l|Kotov fc Vikhlininllioosi ) from XMM-Newton observations 
of clusters at 0.4 < z < 0.7. 

Next, the Yx — M relation was investigated. Yx, the 
product of the gas mass and temperature, has been shown 
to have a low-scat ter scaling relation with cluster mass i n 
simulated clusters (|Kravtsov et al.ll2006l : iPoole et al.ll2007h . 
Yx was calculated for XLSSC 029, and was found to be 
(6.1 ± 1.3) X 1O"M0 keV (using the XMM-Newton kT; the 
Chandra kT gave (7.7 ± 2.2) x IO^^Mq ke V). These values 
are pl otted on the local Yx — M relations of IVikhlinin et al] 
l|2006ll and lArnaud et al.l |20o7) in Fig.d Again, the prop- 
erties of XLSSC 029 are consistent with the local relation. 

Finall y, in Fig. [10] we plot XLSSC 029 on the Lx - M 
relation of lMaugha nl ipOOTl '). This is a sample of 115 clusters 
over the redshift range 0.1 < ^ < 1.3 observed with Chandra. 
For those clusters, the masses were estimated from their Yx 
values. For XLSSC 029, however, we have the advantage that 
the mass was estimated from a full hydrostatic mass analy- 
sis. Once more, the properties of XLSSC 029 are consistent 
with the self-similar evolution of the cluster population. 



6 DISCUSSION 

The deep XMM-Newton and Chandra observations of 
XLSSC 029 have enabled a uniquely detailed study of a 
z > 1 galaxy cluster. This allows the comparison of the prop- 
erties of XLSSC 029 with the expected evolution of the local 
cluster population. The iron abundance in XLSSC 029 is not 
tightly constrained by the data, but the value measured with 




Maughan (2007) data i — • — i 

Maughan (2007) fit 

XLSSC 029 XMM i ' 
XLSSC 029 Chandra i v ' 



1e+15 



M500 (Msol) 



Fig ure 10. XLSS C 029 is plotted in the Lx - M relation 
of jMauehanl [20o3). The two XLSSC 029 points reflect the 
slightly different luminosities measured with the XMM-Newton 
and Chandra data. The central O.lS-Rsoo was excluded from the 
luminosity measurements. 



the central O.lS-Rsoo excluded is consistent with the ob served 
and p redicted evolution in ICM metal abundan ce (|Ettoril 
I2OO5I : iBalestra etHI l2007l : iMaughan et al.ll2008h . The in- 
creased metal abundance when the core regions are included, 
while not statistically significant, is a lso in line with the 
mean trends found bv lMaughan et al] (|2008l ). 

A key result from this analysis is the measurement of 
the temperature profile of XLSSC 029, providing the foun- 
dation for the hydrostatic mass analysis of the system. The 
accuracy of the XMM-Newton temperature profile is sup- 
ported by the agreement of the cruder Chandra profile. The 
modeling of the temperature profile is a key step in the mass 
analysis, and we note that while the model used provides 
a good description of the radial temperature distribution 
in XLSSC 029, other models could also be used to fit the 
data. This is not a concern, as we do not require a phys- 
ical interpretation of the individual model parameters; the 
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model is simply used to give the value and gradient of the 3D 
temperature. Importantly we do not extrapolate our mass 
estimates beyond the range of the data, so the results are 
not strongly dependent on the choice of temperature profile 
model. That said, the model values from beyond the pro- 
jected radial range of the data do have some effect on the 
fit, because gas from all external radii is projected along the 
line of sight in the fitting process. This effect is negligible 
because the low density of the gas at large radii minimises its 
contribution to the projected temperature at the projected 
radii of interest. 

As discussed in SjJlthe analysis of the X-ray images sug- 
gests that XLSSC 029 is a relatively relaxed cluster, par- 
ticularly compared to the high-redshift cluster population. 
The X-ray mass estimate requires that the ICM is in hydro- 
static equilibrium with the gravitational potential. While 
the X-ray morphology suggests that this is the case, there 
are examples of clusters with relaxed X-ray morphologies 
that show evidence for being out of hydrostatic equilibr ium 
(e.g. CLJ1226.9-I-3332 at z=0.89; iMaughan et al.|]2007^ . In 
the case of CLJ1226. 9-1-3332, a temperature map showed 
asymmetric structure coincident with a subclump of galaxies 
detected in the optical images. The X-ray data are insuffi- 
cient for a spectral map of XLSSC 029 (to achieve a similar 
quality map to that of CL J 1226.9-1-3332 would require al- 
most ten times the current XMM-Newton exposure), but 
the available optical data do not indicate substructure in 
the galaxy distribution. 

Many relaxed clusters in the local universe exhibit cool 
cores due to the short rad iative cooling times in the central 
regions fe.o. lFabianlll994l '). The fraction of cool core clusters 
has been found to be significan tly lower at redshifts higher 
than 0.5 (|Vikhlinin et al.ll2007l ). This is believed to be due 
to the higher incidence of cluster mergers at high redshift, 
with mergers disrupting the cooling process. The tempera- 
ture profile of XLSSC 029 gives no indication of any cool 
core in the system. This could indicate that the system has 
recently undergone a merger event, or it could simply be 
because the gas in the cluster core has not had long enough 
to cool significantly. In order to investigate the latter possi- 
bility, the radiative cooling time of the ICM was calculated 
using the measured gas density and temperature profile ac- 
cording to 



8.5 X 10^° yr 



T 



1/2 



,10-3 cm-^ 7 VlO»K' 

l|Sarazin|[l98^ ). and the resulting profile is plotted in Fig. 
111! The profile shows that the radiative cooling time of the 
ICM in XLSSC 029 is longer than the Hubble time (at z = 
1.05) outside of the central ~ 30 kpc, so even without recent 
merger activity, there may not have been sufficient time for 
a significant cool core to develop in this system. 

The gas mass fractions measured for XLSSC 029 at 
-R500 and -R2500 are in good agreement with t he co rrespond- 
ing values measured for the I Vikhlinin et"aLl l|2006l ) clusters. 
This is consistent with no evolution in /gas in our assumed 
ACDM cosmology. With the reliable masses measured in 
these deep observations, XLSSC 029 can provide a useful 
high-redshift datapoi nt for cosmologic al studies that make 
use of /gas evolution (|Allen et allbOO^ ). 

Perhaps the most interesting results from the study of 
XLSSC 029 are obtained by its use as a probe of the evolu- 



Hubbletimeatz=1.05 - 



R (Mpc) 

Figure 11. Cooling time profile of XLSSC 029. The shaded re- 
gion shows the la uncertainties, and the dashed line indicates the 
Hubble time at the redshift of XLSSC 029. 



tion of the mass-observable scaling relations. In all of the re- 
lations tested, the observable properties of XLSSC 029 were 
consistent with the self-similar evolution of the scaling rela- 
tions. In particular, the agreement of XLSSC 029 with the 
Yx — M relation provides further support for the use of Yx 
as a mass proxy for even the most d istant clusters (see also 
iKravtsov erani2006l : lMaughanll2007l '). Furthermore, XLSSC 
029 w as found to lie on the Lx — M relation of iMaughanI 
|2003). This provides additional support for their conclusion 
that Lx can be used as an effective mass proxy, as the mass 
of XLSSC 029 was estimated f rom a full hydros tatic analysis 
rather than from Yx as in the lMaughanI (|2007l l sample. 



A similar study was performed on the z = 0.89 clus- 
ter CLJ1226. 9-1-3332, but that cluste r was found to deviat e 
from the self-similar scaling relations (|Maughan et al.ll2007l ) . 
This was deemed likely to be due to merger activity in the 
that cluster. In the case of XLSSC 029, the current data 
show no evidence for merger activity, but are not sufficiently 
deep to detect temperature substructure like that found in 
CLJ1226.9-h3332. The relaxed appearance of XLSSC 029 
and its agreement with the self-similar evolution of mass- 
observable relations form a self-consistent picture. However 
it is possible that biases on our mass estimate due to a de- 
viation from hydrostatic equilibrium and non-standard evo- 
lution could conspire to cancel one-another out. 



7 CONCLUSIONS 

Our hydrostatic X-ray mass analysis of XLSSC 029 has pro- 
vided a unique opportunity to test the evolution of the mass- 
observable scaling relations ai z = 1. While the strength of 
our conclusions are limited by the use of a single object, we 
found no evidence to reject the simple self-similar evolution 
of the scaling relations. This provides support for the use 
of these scaling relations in cosmological studies with large 
sample of distant clusters where the data do not allow such 
detailed mass analyses. 
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